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ABSTRACT 

Aims. We present an exploration of weak lensing by large-scale structure in the hnear regime, using the third-year (T0003) 
CFHTLS Wide data release. Our results place tight constraints on the scaling of the amplitude of the matter power spectrum 
as with the matter density fim. 

Methods. Spanning 57 square degrees to i^g = 24.5 over three independent fields, the unprecedented contiguous area of this 
survey permits high signal-to-noise measurements of two-point shear statistics from 1 arcmin to 4 degrees. Understanding 
systematic errors in our analysis is vital in interpreting the results. We therefore demonstrate the percent-level accuracy of our 
method using STEP simulations, an E/B-mode decomposition of the data, and the star-galaxy cross correlation function. We 
also present a thorough analysis of the galaxy redshift distribution using redshift data from the CFHTLS T0003 Deep fields 
• ^ that probe the same spatial regions as the Wide fields. 

Results. We find as{flm/0-25)'''^* = 0.785 ± 0.043 using the aperture-mass statistic for the full range of angular scales for an 
^ ^ assumed fiat cosmology, in excellent agreement with WMAP3 constraints. The largest physical scale probed by our analysis 
. - - 1 is 85 Mpc, assuming a mean redshift of lenses of 0.5 and a ACDM cosmology. This allows for the first time to constrain 
cosmology using only cosmic shear measurements in the linear regime. Using only angular scales 8 > 85 arcmin, we find 
O"8(f2m/0.25)°in"^ — 0.837 ± 0.084, which agree with the results from our full analysis. Combining our results with data from 
WMAP3, we find = 0.248 ± 0.019 and erg = 0.771 ± 0.029. 
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1. Introduction 

A primary scientific goal of the Canada-France-Hawaii 
Telescope Legacy Survey (CFHTL^]) is the exploration 
of the properties of the dark matter power spectrum 
and its evolution with redshift using weak gravita- 
tional lensing. The weak lensing signal manifests itself 
in a modification of the apparent galaxy ellipticity in- 
duced by the cumulative weak gravitational shear ef- 
fects of large-scale structure (hereafter cosmic shear). 
The statistical properties of the distortion field, as 
a function of angular scale, reflect the properties of 
the Universe and of the dark matter power spectrum 
projected along the line of sight (see reviews from 
Bartelmann &: Schneider 2001: Van Waerbeke fc Mellierj 



2003t lRefregieij|2003t iMunshi et al.ll2006l) . 

The CFHTLS Deep and Wide surveys have been de- 
signed to maximise the scientific reward of the CFHT 
MegaPrime/MegaCam instrument and in particular to 
produce a high-quality cosmic shear survey. The Deep 
and Wide surveys provide image quality, depth and sur- 
vey size optimised for weak lensing studies as well as 
g' , r', i', z') colo urs over the whol e field to get pho- 



tometric redshifts ( Ilbert et al. 20061 ). Both depth and 



field of view have been increased by roughly one or- 
der of magnitude as compared to the first-generation 
of weak lensing surveys, like the Red Cluster Sequence 
fRCS, fe ockstr a et al.l [2002al) and VIRMQ S-Descart 
( Van Waerbeke et al.ll2000ll200ll . [200^ . 12005?) surveys that 
were carried out at CFHT. 

The first CFHTLS cos mic shear results were ana lysed 
bv lSemboloni et all (|2006l ) and lHoekstra et all (|2006f) who 
demonstrated that MegaPrime/MegaCam provides ex- 
cellent quality data for weak lensing studies. Despite the 
optical distortion of the MegaPrime Wide field corrector, 
after correction, no significant B-modes, nor any obvi- 
ous critical systematic residuals that may affect the weak 
lensing signal at the percent level, were found. The shear 
statistics as a function of angular scale were in good agree- 
ment with the theoretical predictio ns of the most popular 
cosmological models, and Semboloni et al.l (|2006l ) showed 
that the amplitude of shear signal convincingly increased 
with depth, as expected from its sensitivity to redshift. 
These early CFHTLS cosi nic shear data were u sed by 
Semboloni et all (|2006l ) and iHoekstra e t al.' ('2006') to de- 
rive constraints on il^-a^ and bv [Schimd et al. (2007) 
to explore some physical models of dark energy. The re- 
sults were consistent with the past CFHT weak lens- 
ing surveys but their precision was still limited by the 
small sky coverage of the early CFHTLS data and by 
the poor knowle dge of the redshift distribution of sources. 



Benjamin et al.l (^2007) (hereafter B07) overcame these lim- 
itations by using the early CFHTLS Wide data to- 
gether with the Red Cluster Sequence survey, VIRMOS- 
Descart and the Garchi ng-Bonn Deep Survey (GaBoDS, 
iHetterscheidt et"aL 2007) weak lensing surveys, and the 



photometric redshifts of the joi nt CFHTLS-VIMO S VLT 
Deep Survey (VVDS) analysis (jllbert et al.ll2006l) . They 
then derived much more reliable and accurate Qrn-o's con- 
straint, lo wering the up p er lim its on trg to be fully consis- 
tent with lSpergel et al.l (|2007l ). 

The early CFHTLS weak lensing analysis, the joint 
B07 and t he recent Cosmic E volution Survey studies 
(COSMOS, iMassev et alll2007bl) explore only small scale 
lensing. Their cosmological interpretation is therefore sen- 
sitive to the non-linear evolution of the dark matter power 
spectrum and several other physical and systematic effects 
that primarily contaminate the lensing signal at small 
scales. The most serious are the hi gh contribution of non- 
Gaussianity to the error budget (jSemboloni et al.l 120071 ) 
and the signal contamination on s cales below ~20 arcmin 
by the shear-shap e correlation ( Hirata fc SeliakI 2004; 
Hevmans et al.l 200 6b') a nd by the intrinsic ellipticity 
corre lation ( King fc Schneider 2002; Hevmans &: Heavena 
20031) . 

The CFHTLS Wide survey has been designed to probe 
angular scales up to 8 degrees (the largest scale explored 
by all Wide fields). The exploration of angular scales 
beyond one degree is technically challenging due to the 
decreasing amplitude of the lensing signal. Systematics 
in this unexplored territory are also still poorly under- 
stood or unknown. However, they depend on the large- 
scale accuracy and stability of field-to-field astrometric, 
photometric and Point Spread Function (PSF) calibra- 
tions, and thus there is a need for data homogeneity when 
analysing galaxy pairs separated by more than the one- 
degree MegaCam field of view. 

The third release CFHTLS T0003 circumvents these 
issues: The T0003 Wide data explore angular scales up 
to 8 degrees, that is more than one order of magnitude 
larger than the largest non-linear angular scales. It covers 
a total field of view slightly smaller than B07, but with the 
great advantage of forming a single homogeneous sample 
and of being easily calib rated using t he CF HTLS- VVDS 
photometric redshifts of Ilbert et al. ( 2006f ) that are also 
derived from the T0003 release. 

This work presents a weak lensing analysis of the 
CFHTLS T0003 i'-band Wide survey. It extends the pre- 
vious analysis of the CFHTLS Wide to angular scales 
up to 230 arc minutes (about 85 Mpc, assuming — 
0.27 and h — 0.72, a flat Universe and a mean lens 
redshift of 0.5). Its sky coverage is 57 square degrees, 
that is nearly two times larger than early CFHTLS data 
and about 35% of the final CFHTLS wide sky coverage. 
Furthermore, it includes a new uncorrelated field, W2, 
providing a better estimate of the field-to-field variance. 
The shear measurement pipeline is calibrated and its per- 
formance is evaluated using simulated ir nages produced 
by the Shear TEsting Programme rSTEP. lHevmans et al 



2006a; Massev et a l. 2007b). The signal error budget in- 
cludes non-Gaussian corrections to the cosmic variance, 



http:/ /www.cfht .hawaii.edu/Science/CFHTLS / 



using the fitting formulae proposed by ISemboloni et al 
(2007. ). The effective redshift distribution of sources is de- 



L. Fu et al.: Very weak lensing in the CFHTLS Wide 



3 



termined from the CFHT LS T0003 Deep s urvey and cal- 
ibrated using the VVDS (jllbert et al.ll2006l ). 

The paper is organized as follows: In Sect. 2 we give a 
description of the data set, including the image stacking 
procedure used in this work. In Sect. 3 we describe the pro- 
duction of weak lensing catalogues. After a brief review of 
the theoretical background, we present the two-point shear 
results, together with the sky curvature correction needed 
at large angular scales in Sect. 4. The redshift distribution 
is discussed in Sect. 5. In Sect. 6 we show the cosmological 
parameter estimates, discuss the constraints from linear 
scales and compare to other data sets. In Sect. 7 we dis- 
cuss the contamination to our weak lensing measurement 
from shear-shape correlations. Finally, we summarise and 
give our conclusions in Sect. 8. 



2. Data Description 

2.1. Overview of the CFHTLS T0003 release 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) is a 5-year project set up jointly by the 
Canadian and French agencies. The Deep and Wide obser- 
vations are all carried out in service mode by the CFHT 
operation staff using the MegaPrime/MegaCam instru- 
ment mounted at the prime focus of the telescope. The 
MegaCam camera is composed of an array of 9 x 4 CCDs 
(2048 X 4612 pixels each). The pixel size at MegaPrime 
focus is 0'.'186, so that MegaCam compris es a compact 
field of view of 1° x 1° (|Boulade et al.ll2003l ). 

Details o n the Deep and Wide fields have been in- 
troduc ed in ISemboloni et al. ( 2006h and Hoekstra et al 



( 2006t ). respectively. After completion the Wl, W2 and 
W3 Wide fields will be composed of 8x9, 7x7, 7x7 dif- 
ferent MegaCam pointing positions, respectiveljU. Each 
centre position is separated by its nearest neighbour fields 
by about one degree. For each field, a sequence of 7x620- 
second i'-band exposures, separated by a small dither, 
is taken. The dithering pattern is encompassed within a 
3' x 4' box. Hence, neighboring pointings overlap in right 
ascension by a minimum of two and a maximum of three 
arc minutes, whereas the overlap in declination is bounded 
between three and four arc minutes. The overlapping re- 
gions are used for the pointing-to-pointing internal astro- 
metric calibration and flux-rcscaling processes. 

The CFHTLS Wide T0003 release is produced from 
all MegaCam CFHTLS images obtained between June 
1st, 2003 and September 5th, 2005, that passed both 
the CFHT and initial Terapix validation processes. 
Each individual raw image has been pre-processed 
(bias/dark/fringe subtractions and flatfielding) , CCD- 
to-CCD flux-rescaled and photom etrically calibrated at 
CFHT using the Elixir pipeline i Magnier &: Cuillandrd 
20041 ). The Elixir products are archived at CADC and then 
transferred to the Terapix data centre for further higher 



level processing and the production of the CFHTLS re- 
leasefl 

In this paper we use the i'-band pointings from the 
T0003 release, but we introduce a more severe image se- 
lection in order to optimise and homogenise the depth and 
the image quality over the whole field used for weak lens- 
ing studies. Each i'-band image must fulfil the following 
second-level criteria before entering into the calibration 
and stacking processes: 

— half-light diameter of the Point Spread Function below 
0'.'9; 

— individual exposure time larger than 500 seconds; 

— at least 4 exposures per stack. 

Ficld-to-field seeing and depth variations are then min- 
imised. Over the 63 i'-band pointings, 57 pass the second 
level weak lensing criteria and have been included in this 
work. They are distributed as follows: 19 pointings on Wl, 
8 on W2, and 30 on W3 (see Fig.[l]). The W2 field is signif- 
icantly less covered than Wl and W3 but is useful in order 
to derive an estimate of the CFHTLS Wide field- to- field 
variance. 



2.2. Image Production 

The input stack images used in this work are produced us- 
ing the same Terapix procedures and software tools as the 
T0003 CFHTLS release. Terapix first generates individual 
weight map images and individual primary catalogues. It 
then proceeds to astrometric calibration, MegaCam field- 
to-field photometric rescaling, image re-centering, image 
resampling, image warping and, finally, image stacking. 
Both a co-added image and its weight map are produced, 
as well as an ASCII DS9 readable mask file and a se- 
ries of quality control meta-data. A description of these 
processing steps and software tools can be found on the 
Terapi x web sitjl as well as on the Terapix release doc- 
ument (jMellier et al.ll2005h . All T0003 configuration files, 
parameter lists and processing command lines are archived 
at CADC. Only i'-band images are considered in this work 
since other filters cover a much smaller field of view, with 
a large scatter in sky coverage and depth between each 
filter. 

All fields are astrom etrically calib r ated and fiux- 
rescaled using SCAMfH (|Bertinl l2005al l2006l) . The as- 
trometric reference catalogue is USNO-Bl, which is suf- 
ficiently accurate for the external astrometric precision 
needed for this work. Internal astrometry and MegaCam 
pointing-to-pointing flux rescaling is done by identifying 
common objects located in each overlap area. The im- 
age re-cente ring, resampli ng and stacking are produced 
by SWarpi ()Bertinll2005bh . using the same configuration 



Details on the T0003 release can be found at 



^ http : //terapix . iap . f r /cplt/oldSite/Descart7] 
summarycf htlswide . html 



http: //terapix. iap. f r/rubrique .php?id_rubrique=208] 


http://terapix.iap.fr/soft 




^ http : //terapix . iap .f r/sof t/ scamp 


® http : //terapix . iap .f r/ soft/ swarp 
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Fig. 1. Sky coverage of the Wl, W2 and W3 fields used 
in this work. Each CCD is drawn as a small rectangle and 
each MegaCam field is a squared mosaic of 36 rectangles. 
The small white holes are regions with missing data. 



and image proce s sing p arameters as those discussed by 
McCracken et al.l (|2003f ). All stacked images have a pixel 
size of 0'.'186. 

The astrometric calibration was performed for each 
pointing individually. For each pointing, only exposures 
located inside a circle of radius of 1.5 degrees were se- 
lected. This circle encompasses all images at the centre 
field position together with all exposures located around, 
at the 8 nearest neighbour centre positions. We experi- 
enced that selecting images located beyond this radius did 
not improve the accuracy and stability of the astromet- 
ric solution, and sometimes would have even degraded it. 
The 8 nearest neighbour fields provide enough common 
stars in overlap regions to stabilise the solutions at the 
boundary of each field. The internal rms error estimates of 
the astrometric calibration derived from the common ob- 
jects of nearest neighbour fields is 0'.'030± O'.'OlO for both 
MegaCam directions, where the uncertainty denotes the 
mean field-to-field scatter. The external rms error is to- 
tally dominated by the USNO-Bl internal error, which is 
0'.'35 in both directions. 

In contrast, each stack does not use nearest neighbour 
images, but only composes together a sequence of expo- 
sures having a small dither with respect to a centre field. 
Each stack is produced by SWarp, using the weighted me- 
dian value of each pixel and a LanczosS interpolation ker- 
nel. All output images have 19354 x 19354 pixels of 0'.'186, 
with North-East orientation along the X and Y pixel co- 
ordinates. For all images we use a tangent projection and 
the Equatorial J2000.0 astrometric coordinate system. 

The reference photometric zero-point has been de- 
rived by CFHT using standard Landoldt calibration fields 
(Landoldt 1992), but all catalogues produced prior to 
weak lensing analysis have a default zero point magni- 
tude set to 30.0. The magnitude system is instrumen- 
tal AB. An inspection of stellar colour-colour diagrams 
of each field observed in 5 bands shows that the field- 
to-fleld scatter in the overlapping regions is 0.03 mag. 
Comparison of SDSS and CFHTLS common stars shows 
that the z'-band photometry agrees within 0.01-0.02 mag- 
nitude rmjll- However, only 10 W3 and 2 Wl fields have 
common objects with SDSS, so similar external quality 
assessments cannot be done for all Wide pointings. 

The mask files produced at image processing consist 
on a set of polygons defined for each pointing in WCS 
coordinates. They mask the periphery of each MegaCam 
field of view and all halos and saturated spikes produced 
by bright stars. In order to avoid contamination by halos 
or diffusion from very luminous objects, all bright stars 
located up to a radius of 45 arc minutes from the centre 
position are automatically masked. The size of polygons is 
scaled to their apparent magnitude provided in the USNO- 
Bl catalogue. 

The masks are then tuned by adding or modifying 
polygons from a visual inspection of each stacked image. 
This step is necessary to clean all images from non-stellar 



^ http : //terapix . lap .f r/eirticle ■php?id_article=593| 
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contamination or stellar defects that were missed by the 
automatic masking process. This includes big halos pro- 
duced by extremely bright stars, nearby galaxies or any 
features that may produce a diffuse light component with 
sufficiently steep gradient to contaminate the second mo- 
ments of a galaxy's surface brightness profile which is used 
to derive its ellipticity. Regions with low signal-to-noise 
ratio are also masked. In particular, the imprints of gaps 
between CCDs as well as the boundary of each field are 
discarded and masked systematically. They are revealed 
by low-noise strips with a typical rectangular shape that 
draws the border of each detector. Finally, meteorite, as- 
teroid and aeroplane tracks that may still remain in the 
stacks are masked as well. The size of each polygon is gen- 
erally significantly larger than the visual size of the defect 
it masks. Using this conservative masking process, the fi- 
nal effective sky coverage of the 57 selected Wide fields 
drops to 34.2 deg^, roughly 60% of the total field. 

3. Production of weak lensing catalogues 

Our shear measurement pipeline was optimised and cal- 
ibrated using the STE Pl and STEP2 s i mulatio ns from 
Hevmans et al. ( 2006a ) and Massev et al. ( 2007bl ). See ap- 
pendix A for a description of both our pipeline and the 
STEP simulations. Table ETT] lists the optimised parame- 
ter values of our pipeline. 

3.1. Object selection 

The lensing catalog ue is generated by the IMCAT software 



= 0.1 



([Kaiser et al.lll995f l. The size of each object is defined by 



the aperture radius parameter given by the IMCAT peak 
finding algorithm. The significance detection threshold, as 
defined by the IMCAT parameter ly, is set to = 8 {i.e. 
above the rms noise). This value was set according to the 
STEP tests in order to maximise the number of objects 
detected while still keeping the bias on the shear compo- 
nents negligible. The catalogue is then filtered to remove 
objects with radius smaller than the seeing or larger than 
6.75 pixels (about 1.3 arc seconds). Pairs with angular sep- 
aration smaller than 10 pixels (1'.'86) are also discarded in 
order to avoid contamination from overlapping isophotes. 

The magnitude of each object is derived by computing 
its fiux within an aperture radius of 3 x r^. Only objects 
with IMCAT magnitude 21.5 < i'^g < 24.5 are kept into 
the final analysis catalogue. Beyond this limit, the sample 
completeness drops significantly below 50%, most objects 
are too noisy and their shapes are no longer reliable for 
the precision needed for weak lensing studies. The final 
catalogue based on the T0003 release of CFHTLS Wl, 
W2 and W3 fields contains roughly two million galaxies. 
Due to the different weighting applied during the sam- 
ple selection, the effective number of galaxies used for the 
weak lensing analysis is 1.7 million, spread over the effec- 
tive area of 34.2 deg^. It corresponds to a galaxy number 
density of 13.3 gal/arcmin^. The shapes of these galaxies 
are quantified by measuring their ellipticities. 
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Fig. 2. The pattern of the PSF anisotropy in an example 
pointing W3-t-2-h0-CFHTLS_WJ.143023-f 543031. Ticks 
represent the observed ellipticities at stellar locations. On 
top of the figure a 10% ellipticity modulus is shown for 
comparison. 



3.2. PSF correction 

The ellipticities of galaxies are corrected from the 
PSF produced by telescope, detector, optical and 
atmospheric effects , using the Ka i ser, Squires and 
Broa d hurst method ( Kaiser et al ] ll995l : iLuppino fc Kaisen 



19971: IHoekstra et all I1998I) . hereafter KSB+. Our im- 
plementation of KSB-I- is based on the o ne used in 
Van Waerbeke et al. I (I2OOOL I2OOIL I2OOI I2OO5I . referenced 
as 'LV in Heymans et al. 2006a). We calibrated it and 
modified its input parameters after a new series of opti- 
misations made with the STEP simulations. The results 
are presented in Appendix A. The version used in this 
work recovers shear with an underestimation of only 1% 
to 3% on the simulated images. 

The PSF is measured at the stellar positions. After 
identifying the stars in the images and assuming the PSF 
changes smoothly across the field, the KSB-quantities 
known at the stellar positions can be estimated at the 
galaxies positions by using a polynomial fit. The typical 
pattern of the PSF anisotropy across one 1 deg^ field shows 
a significant variation across the whole camera (Fig. [2]) 
which suggests the need to perform the fit in each CCD 
separately. Each CCD covers 7x4 arcmin^ and contains 
an average of 43 stars, which allows an accurate mapping 
of the PSF with a second order polynomial function. 

A weight, w, is assigned to the ellipticity components 
of each galaxy and us ed in the shear measur ement (see 
Eq. [TU]). We use the iHoekstra et~al] (|2002bf ) weighting 
scheme 

p2 

(1) 
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Fig. 3. The average galaxy weight (with arbitrary normal- 
isation) as a function of i'j^g in the range of [21.5;24.5]. 



where gp is the error on the e llipticity measurement de- 



fined in Hoekstra et al.l and Pj is a shear polaris- 

ability (jLuppino fc Kaiserlll997l ) . The weight also depends 
on the intrinsic ellipticity dispersion, Ce, which is derived 
from the mean intrinsic ellipticity dispersion of the whole 
galaxy catalogue. We find a value of ae = 0.42. 

The shape of the weighting as function of the magni- 
tude is shown in Fig. [3l It decreases for fainter magnitudes 
since the error on the ellipticity increases when the signal- 
to-noise ratio decreases. 

4. Two-point cosmic shear statistics 

4.1. Theoretical background 

The cosmic shear power spectrum is identical to the lens- 
ing convergence power spectrum, P^, which is a projection 
of the dark matter po wer spectrum, Pa, along the lin e of 
sight (see for example iBartelmann fc Schneideiil200lh : 



-ni 



dx 



Xli. 



Mx' ~ X) 



fKix') 



fKixY 

2 



(2) 



where x is the comoving distance along the light ray and 
Xiim is the limiting comoving distance of the survey; Jk (x) 
is the comoving angular diameter distance; n{x) is the red- 
shift distribution of the sources and £ is the modulus of a 
two-dimensional wave vector perpendicular to the line of 
sight. Equation ([2]) shows that the cosmological informa- 
tion contained in the lensing power spectrum is degenerate 
with the redshift of the sources. 

The convergence power spectrum can be derived from 
the two-point shear correlation functions. In particular, 
the ^± correlation functions relate to the power spectrum 



according to 

^±{0) = Ctt{0)±^xx{0)^^ deiP,{e)3oAm , (3) 

where ^tt(^') and x (^) are the tangential and rotated el- 
lipticity correlation functions (given in EqllOp , 9 is the an- 
gular separation between galaxy pairs, and J0.4 are Bessel 
functions of the first kind. 

Other two-point functions of the shear field may be 
derived from such as the top-hat filtered variance of 
the shear a nd the variance of the a perture- mass, in circular 
apertures ( Schneider et al. 2002bl ). Respectively, 

(|7P>E3W = 



dt9i? 



(4) 



and 



29 



/o 29^ 
The filter functions 5*. 



e+(^?) ± T_ 



(5) 



Schneider et all (|2002bD . 
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xy/A-x^{6 - x^) - 8(3 - x^) arcsin(x/2^ 



(6) 



T+{x) = 



6(2- 15a;2) 



■ arcsin(a;/2) 



r;%/4~ 



IOOtt 



(120 + 2320x2 - 754x'* -I- 132a;^ - 9a;^) H(2 - x) 



T-ix) 



192 
35^" 



x^ I 1 



X 

T 



2\ 



H(2 - x) , 



(7) 



where H denotes the Heaviside step function. 

All second-order statistics are different filtered versions 
of the convergence power spectrum. Therefore they probe 
different properties of the same power spectrum. 

The cosmological shear field is (to first order) curl-free 
and is called an E-type field. It is useful to decompose the 
observed shear signal into E (non-rotational) and B (ro- 
tational) components. A detection of non-zero B-modes 
indicates a non-gravitational contribution to the shear 
field, which reveal s a likely system a tic co ntam ination to 



the le nsing signal. ICrittenden et al.l (|2002f ) and iPen et al 



(|2002l ) derived an analogous decomposition for the shear 
correlations, which is also used in this work: 



&,B(e) 



(8) 



where t he definition of ^' is also given in Schneider et al.l 
()2002bl ) 



(9) 
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Both S- and ^' have infinite support, which imphes the 
E/B decomposition of the shear correlation function and 
of the top-hat shear variance are non-local. They can be 
computed from data, up to an offset value which depends 
on the largest angular separation i9max- This offset is a 
constant for ^e,b and a function of 6 for {\j\'^)e.b{S)- In 
contrast, the aperture-mass variance decomposition is lo- 
cal, providing an unambiguous decomposition. In practice, 
however, a lower limit on the angular separation of galaxy 
pairs imposed by contamination of overlapping isophotes 
may bias its amplitude. The lack of galaxy pairs closer 
than around 3" causes an underestimation of the aperture- 
mass dispersion. This bias is small, of order 5% for 9 — 1' 
and smal ler than one percent on scales larger than 2 arc 



minutes ( Kilbinger et al. 20061 ) 



4.2. Sky curvature correction at large angular scales 

The shear correlations are computed as follows: 
J2wiW.jet{:x.i) ■ et(xj) ^ 



J2 ^i'^l 

J^WjWjexi^x^i) ■ ex(xj) 



(10) 



where 9 the separation of pairs. The elliptici- 

ties are locally decomposed in each pair frame in a tangen- 
tial and a cross-component. The tangential component is 
computed orthogonal to the line connecting each galaxy 
pair. The cross-component is derived at a tt/4 angle to 
the connecting line. Using Eqs. ([4])- ([8]) we estimate the E 
and B modes of the top-hat variance, of the aperture-mass 
variance and of the shear correlations. 

We correlate galaxies which are up to more than seven 
degrees apart. At such large angles the curvature of the 
sky is no longer negligible. To avoid a potential bias due to 
projections we calculate distances and angles in spherical 
co-ordinates as follows. 

The distance d between two objects at right ascen- 
sion and declination {ai,Si),i = 1,2, computed along the 
great-circle connecting the 2 objects is given by 

cosd — cos(ai — 02) cos^i cos 52 + sinSi 811162. (H) 

In order to decompose the ellipticity in tangential and 
cross components, we need to consider the angle between 
the great circle defined by the two galaxies and a paral- 
lel, since each pointing is aligned with lines of constant 
declination. This is the course angle, given by 



smfai — 02) cos (52 
tan ip — - 

cosoi sin02 — sinoi COSO2 cos(q!i — a2) 



(12) 



Not only the distances but also the galaxy ellipticities 
are affected by the sky curvature. The ellipticity modu- 
lus of a galaxy remains virtually unchanged, curvature on 
scales of a few arc seconds is negligible. However, since 
ellipticity is also characterized by its orientation, one has 
to be careful when correlating ellipticities at large angular 



distances (|Castro et al.l 120051 ) . In our case, the ellipticity 
components ei and 62 are measured in the local Cartesian 
co-ordinate frame given by the X- and Y-axes of the cor- 
responding individual CFHTLS Wide pointing. Since each 
pointing is projected using its own tangent point and de- 
fines its own co-ordinate frame, we correlate ellipticities 
of galaxies from different pointings by using their com- 
ponents measured in the respective local co-ordinate sys- 
tems. In doing that, the sky curvature is neglected over 
the scale of a single pointing but it is taken into account 
between pointings. The effect in a single pointing corre- 
sponds to a small ellipticity rotation for galaxies which 
are at a finite distance from the pointing centre, at most 
43 arc minutes. 

We compared the shear statistics computed using 
spherical co-ordinates and using the following simple pro- 
jection: Cartesian co-ordinates {X, Y) of a galaxy with 
right ascension a and declination S are defined by X = 
a cos (5c and Y = S, where 6c is the declination of the field 
centre. The relative error is on the order of a couple of per- 
cent on average. On larger scales, where the shear signal 
is small, this relative error can be much higher. Therefore, 
throughout this paper we take the sky curvature into ac- 
count by calculating the shear statistics in spherical co- 
ordinates. 

4.3. Results 

The shear correlation functions ^± are computed in narrow 
bins. We use angular separations in the range between a 
conservative lower limit of 3 arc seconds and a maximum 
separation i?max where the number of pairs per bin be- 
comes very small. For each field, the number of pairs per 
bin shows a similar "top-hat" behavior: a very steep in- 
crease from -d = followed by a roughly constant value 
up to a i9max where it starts a very steep decrease to zero. 
This separation iJ^ax is 400' for Wl, 240' for W2, and 462' 
for W3. 

From the two-point shear correlation functions we cal- 
culate the shear top-hat variance and the aperture-mass 
dispersion up to a radius of 0max which is half of the largest 
separation z?max, according to Eq. ([4]) and ([5|). 

The missing information for ^e,b(^) and (|7P)e,b(^) 
on scales larger than 20niax is accounted for by adding 
theoretical predictions of these off-sets to the data using 
a fiducial cosmological model. Alternatively, we may set 
the B-modes of the shear correlation function and top- 
hat variance to zero on the angular scales where we mea- 
sure zero aperture-mass dispersion B-modes. We checked 
that both methods produce very similar and small off- 
set values and thus this procedure does not bias the cos- 
mological interpretation towards the fiducial model used. 
Furthermore, our cosmological estimates are made using 
the aperture-mass dispersion and are free of this small 
arbitrariness. 

The three statistics are plotted in Fig. |4] and the corre- 
sponding values and errors are provided in Table IB.1IIB.2I 
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Fig. 4. Two-point statistics from the combined 57 point- 
ings. The error bars of the E-mode include statistical noise 
added in quadrature to the non-Gaussian cosmic variance. 
Only statistical uncertainty contributes to the error bud- 
get for the B-mode. Red filled points show the E-mode, 
black open points the B-mode. The enlargements in each 
panel show the signal in the angular range 35'-230'. 
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Fig. 5. The top-hat E-mode shear signals of Wl up to 
200', of W2 up to 120' and of W3 up to 230' are shown. 
The error bars includes statistical noise and cosmic vari- 
ance for each individual field. 



in Appendix B. It is worth noting that this is the first 
time that a cosmic shear signal has been measured down 
to i'j^g = 24.5, beyond scales of one degree . Notice also 
that the independent measurements of the shear statis- 
tics made in the three individual fields Wl, W2, W3 are 
statistically consistent at all scales. This is illustrated by 
Fig. [5l where the three measurements of top-hat disper- 
sion are shown. 

In Fig. [4] the error bars of the E-modes include sta- 
tistical noise and cosmic variance calibrated for non- 
Gaussianity, while the error for the B-modes only includes 
statistical uncertainty. We find a clear E-mode signal and 
a B-mode which is consistent with zero throughout the ex- 
plored range of angular scales, except between 50 and 130 
arc minutes where there is a small but significant feature 
in all three second-order functions. This bump of the B- 
mode peaks at about 60-80 arc-minutes which are the side 
and diagonal sizes of a Megacam field. We therefore guess 
it is due to a correlation in PSF residuals on the scale of 
the camera. In Sect. |6]we show that our cosmological re- 
sults are not biased by this level of residual systematics 
on this range of angular scales. 

On very large scales (120'-230') we find a very small 
B-mode, much smaller than both the E-mode amplitude 
and cosmic variance, but which is not always within la 
of a zero detection. Notice that the errors on the B-mode 
shown in Fig. [4] are theoretical (statistical) and not es- 
timated from the data, which would include systematics 
(for example error contributions may arise from the in- 
complete PSF correction). Moreover, the signal-to-noise 
with the present CFHTLS Wide data is so high, even for 
B-modes, that subtle effects may dominate the very small 
Poissonian error, particularly on large scales where there 
are a significant number of galaxy pairs. 
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Fig. 6. Left panel: Binned scatter plot of the shear estimates (one component) using the two pipeUnes. Dark colours 
show highest density of points. The bin size is 0.05 in ei. Right panel: The aperture-mass variance from Wl, W2 
and W3 measured with the two pipelines up to scale 210 arc minutes, using only objects which are detected in both 
pipelines. For clarity of the comparison, the error bars only show the statistical errors, but the cosmological analysis 
of this work includes the whole error budget (see text). These error bars are larger than the one of Fig. 0] because 
number of common objects are smaller the full catalogue. Note that the large negative B-mode on small scales is not 
present in the full catalogue, see Fig. [4l 



The field-to-field variation of the B-modes is a possible 
way to assess these effects on the error buget. We tried to 
measure this by splitting the 3 Wide fields into 1 1 blocks of 
2x2 deg^ each, which allows to calculate the B-modes on 
scales up to 60 arcmin in each block. We obtained B-modes 
with amplitude very similar to Fig. |4] but the field-to-field 
scatter is larger than the plotted error bars and reaches 
a factor of 2 at 60'. This is an interesting indication that 
we are likely underestimating the error on B-modes, even 
though it is not a precise measurement due to the small 
number of independant fields. A thorough analysis of this 
noise contribution needs many more field and is left to a 
future analysis of the CFHTLS four year data. 

4.4. Cross-check and control of system atics 

We cross-checked the shear measurement by using an inde- 
pendent analysis on the same data sets. This analysis was 
done with another version of KSB-I- that has been tested 
with the STEPl-l-2 simulations ('HH' in iHevmans et al 



l2006at iMassev et al. 2007bf ). Hereafter, we refer to our 
analysis as 'Pipeline F and to the 'HH' results as 'Pipeline 
H'. 

The left panel of Fig. [6] shows the shear estimated for 
each galaxy by each of the pipelines. The results are in 
good agreement for ellipticity values per component be- 
tween —0.6 and 0.6. For ellipticities outside this range the 
dispersion between the pipelines is larger and a trend for 
an underestimation of the shear from Pipeline I with re- 
spect to Pipeline II can be seen. Note however that the 



pipelines are not optimised for large ellipticities, since 
the STEP simulation galaxies have ellipticities that are 
smaller than 0.1. 

We then compare the two-point functions using the 
aperture-mass variance. We choose this statistic because 
angular scales are less correlated than for the top-hat dis- 
persion. Moreover, it does not have any ambiguity re- 
lated to a non-local E/B decomposition. The values of 
Map are calculated from the two pipelines using only ob- 
jects detected by both pipelines. Because the pipelines 
have different selection criteria the common objects are 
only two-thirds of the whole sample. Each object is as- 
signed a weight which is the product of its weights in each 
of the two pipelines. The largest radius explored in the 
comparison is 210 arc minutes. As can be seen in Fig. [S] 
(right panel), the E- and B-modes of the two pipelines 
are remarkably similar. The differences are within the la 
errors on all angular scales. The small B-mode bump ap- 
pears in both results at 60-80 arc minutes, as in Fig. H) It 
also drops to nearly zero at all scales beyond 120 arcmin 
for both pipelines. The bias between the two pipelines at 
large ellipticities, seen in the left panel of Fig. [6] is not vis- 
ible here. The reason is that the large ellipticity galaxies 
represent less than 4% of the sample. Furthermore these 
galaxies are typically downweighted; large ellipticities are 
difficult to measure, resulting of a larger error on the ellip- 
ticity measurement, and the shear polarisability increases 
with ellipticities. They have consequently a lower weight 
according to Eq. |T]). 
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These results are not expected to be identical to the 
aperture mass dispersion computed with the whole sam- 
ple, shown in Fig. 31 because the number of objects in 
the two samples is different. They are however similar, 
except for the large B-modes on scales smaller than 2', 
which are detected by both pipelines on the smaller sam- 
ple. Since both analyses use KSB, these B-modes may be 
due to similar residuals of the PSF correction, but we can- 
not rule out an intrinsic B-mode contribution. Whatever 
the origin we only use angular scales larger than 2' for the 
cosmological parameter constraints (see also Sect. 14. ip in 
order to avoid any contamination. 

The most common and problematic source of con- 
tamination of the lensing signal is the imperfect PSF 
anisotropy correction. The angular dependence of any PSF 
systematic residual may be checked by computing the cor- 
relation between the correc ted galaxy a. n d un corrected 
stellar ellipticities. Following i Bacon et al. I (|2003 l) we nor- 
malise this quantity by the star-star uncorrected ellipticity 
correlation in order to assess its amplitude 



(e*(x)7(x + g))^ 
{e*{x)e*ix + 9)) 



(13) 



where the symbol -k indicates a stellar quantity. Figure [7] 
shows this cross-correlation compared to the shear signal 
up to 230'. Overall, the amplitude is at least one order 
of magnitude smaller than the signal. This demonstrates 
that the PSF correction is under control in our shear anal- 
ysis up to an angular scale of nearly 4 degrees. 

5. Source redshift distribution 

The calibration of the source redshift distribution in the 
the CFHTLS Wide fields cannot be calculated from the 
Wide photometric data since only a few fields have al- 
ready been observed in 5 bands. However, the CFHTLS 
Deep fields overlap, or are located very close to, the Wide 
fields. One can therefore use the p hotometric re dshifts de- 
rived for the CFHTLS Deep data (jllbert et allEoO^ as a 
representative sample of the Wide galaxy population, in 
particular for Wl that covers the Dl field. 

5.1. Building of the parent Deep n{z) liistogram 

The lllbert et al.l (j2006t ) catalogue samples photomet- 
ric redshifts of more than 500,000 objects in the four 
CFHTLS Deep fieldfH, with an limiting magnitude 
much fainter than that of the Wide survey, covering the 
range < z < 6. It has been calibrated with spectro- 
scopic redshifts obtained by the VVDS Surve y in the 



CFHTLS Deep Dl field (|Le Fevre et al.ll2005f ). In this 



photometric redshift catalogue 318,776 galaxies have a 
magnitude matching the range used in the Wide survey, 
i.e. 21.5 < i'j^g < 24.5. This sub-sample is used to build 
up our redshift distribution. 



For each object in lllbert et al.l (|2006l ). the released pho- 
tometric redshift catalogue provides the maximum likeli- 
hood redshift Zp and error estimates such as the left and 
right la error. In order to estimate the redshift distribu- 
tion we build a normalised Gaussian probability distribu- 
tion for each galaxy, with mean Zp and dispersion given 
by the mean of the left and right error. We then draw a 
redshift z randomly and repeat the procedure 1000 times. 
The variance of these 1000 randomizations is considered 
into the final error budget. 
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Fig. 7. The cross-correlation function ^sys (Eq- fT3)l be- 
tween galaxies and stars is shown as a function of angu- 
lar scale up to 230' (black filled). The amplitude of the 
cross-correlation is always at least one order of magnitude 
smaller than the shear signal (red open). 



5.2. Rescaling to the Wide population 

To take into account the different selection functions be- 
tween the Deep parent sample and the Wide catalogue 
used in this work, each galaxy is weighted according to 
the ratio of the Wide to Deep galaxy number density, see 
Fig. [S] In addition, we include the weak lensing weight 
(Fig. [3]) to match the redshift distribution to the weighted 
galaxy population selected for weak lensing. The redshift 
distribution is built up with all photometric redshifts in 
the range < z < 4. 

5.3. Error budget 

The errors on the histogram have several contributions. 
First, the uncertainty in the photometric redshifts is es- 
timated from the variance of the 1000 randomizations 
from the CFHTLS Wide redshift histogram constructed 
in Sect. 5.1. Second, Poisson noise, crp is added as ^/n, 
where n is the number of galaxies per redshift bin. Third, 
we need to include sample variance, CTsv, since we estimate 
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Fig. 8. Magnitude distributions for the Deep (solid line) 
and Wide (dotted). We use the ratio Wide/Deep for the 
rescaling of our redshift distribution. The Wide effective 
number density takes into account all weak lensing selec- 
tion criteria and has been multiplied by their correspond- 
ing weights (Fig. [3]). 
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Fig. 9. Final normalised redshift distribution. Galaxies 
are selected in the range [0;4], and the best-fit is given 
for function given in Eq. ([M]) . Note that the fit is only 
performed in the interval [0;2.5]. 



the redshift distribution from a reference catalogue. The 
sample variance is given as a functio n of Poisson noise and 



redshi ft for various survey areas in iVan Waerbeke et al 



We use the dsv/cp ratio of a one square degree 
survey, corrected for our bin size. We further rescale it ac- 
cording to the weak-lensing selection function, since this 
reduces the total number of galaxies, on which the ratio 
depends, /ap cx nj yjn. Note that we do not divide 

the ratio by \f\ to account for having four independent 
Deep fields, since the Poisson error is calculated for the 
sum of the four fields. 

With the large number of galaxies in our sample and 
the high accuracy reached by the photometric redshifts at 
z 1, the sample variance is the dominating contribution 
to the error budget. Poisson noise and redshift uncertain- 
ties only contribute ~ 5% at z = 1 but become dominant 
for z > 3 where the number of galaxies is very low. As a 
cross-check, we have calculated the field-to-field variance 
of the four Deep photometric redshift catalogues. The re- 
sult is consistent wi t h the sample variance obtained by 
Van Waerbeke et al. using numerical simulations. 



ergy distributions of galaxies with z < 0.2 and z > 1.5. 
The recent analysis of the sp atial correlation of popula- 
tions in different photo-z bins (jVan Waerbeke et al.ll2007 . 
in prep.) confirms that more than 50% of galaxies in the 
peak are most probably at redshift z ;$ 0.4. 

We do not have a reliable estimate of the histogram 
bin-to-bin correlation. Indeed, the off-diagonal sample 
variance was not calculated in the n umerical simulation 
analysis of Van Waerbeke et al. ( 20061 ). and a field- to-field 
estimate using the four Deep fields is too noisy to be of 
practical and reliable use. Thus, in order not to propagate 
systematics present in the histogram into the cosmologi- 
cal constraints it is preferable to use a fitting function to 
the redshift distribution in the cosmological parameters 
estimation. For this we consider all galaxies in the range 
< z < 2.5 and fit the redshift distribution with the 
following function. 



n{z) = A 



z" + z°-^ 



A = 



dz ) \ (14) 



5.4. Fitting n{z) of the Wide weal<-lensing sample 

A histogram of the sources redshifts is shown in Fig. [51 
where the error bars include redshift uncertainty, Poisson 
noise and sample variance. Although sample variance is 
taken into account, the histogram shows a significant 
bump at redshift z ~ 3. We cannot exclude the possi- 
bility that this small peak might be partly a real feature 
resulting from the joint spectroscopic, photometric and 
weak lensing selection functions of our galaxy sample. It 
is however more likely to be an artifact due to system- 
atic photometric redshift misidentifications arising from 
degeneracies that exist between the optical spectral en- 



Thc normalization A is determined by integrating until 
^max — 6, the upper limit of the photometric redshift cat- 
alogue. 

This function provides a better fit to the main peak 
and the tail of the distribution as compared to the power- 
law function used in B07. The distribution shown on Fig. [9] 
corresponds to the best-fit parameters listed in Table [TJ 
As expected, the peak at z '--^ 3 is no longer present. It 
is worth mentioning that, although the histogram shows a 
significant fluctuation with respect to the best-fit model at 
redshift z ^ 3, the mean redshift derived from the best-fit 
distribution is within 1% of the mean value of the his- 
togram. 
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Table 1. Results of the fit to the redshift distribution 
< z < 2.5, using Eq. ([M]) . The la error bars of three 
parameters are shown as well, (z) is the mean, 2:,„ the 
median redshift. 



a 


0.612 ± 0.043 


A 


1.555 


b 


8.125 ±0.871 




0.949 


c 


0.620 ± 0.065 




0.826 



6. Cosmology with CFHTLS Wide 

6.1. Shear covariance 

The covariance matrices for the shear two-point corre- 
lation functions are ca lculated using the expressions of 
ISchneider et al. ( 2002a ). valid for a Gaussian shear field. 
They consist of a statistical noise term, a cosmic variance 
term and a mixed term. To account for possible resid- 
ual systematics in the shear signal, we add the measured 
B-mode at a given angular scale quadratically to the cor- 
responding diagonal element of the covariance. 

The first three terms are calculated using a Monte 
Carlo method applied to the measured galaxy positions 
and their weight similar to the bootstrapping defined 
in Sect. 5. In that way the survey geometry, bound- 
ary effects and the non-unif orm, discrete galaxy distribu - 
tion are taken into account (jKilbinger fc Schneideiil2004l ) . 
Furthermore, this method allows to compute a statisti- 
cal noise that not only includes the shape noise of the 
two-point functions estimators but also takes into account 
Poisson or shot noise. 

The non-Gaussianity of the shear field on small 
scales is considered by applying a correction to the 
cosmic variance term using the calibration formula of 



Semboloni et al.l (|2007l ). The parameters for the model 
shear correlation function are Om = 0.27, f^A = 0.73, h = 
0.7, Qb = 0.044 0-8 = 0.8 and = 1.0, using the 



Smith et al. I (|2003l . hereafter SOS) non-linear prescription. 



The redshift distribution is the best-fit of the n{z) data 
(see Sect. 5). For the non-Gaussian calibration a mean 
redshift of 0.95 was assumed. 

The top-hat variance, the aperture-mass statistic and 
the E-/B-correlation functions are functions of both ^+ 
and (Eqs. HUH]). Therefore, their covariance matrices 
depend on the full covariance of the combined data vector 
(^+,^_). However, we use only C++, the covariance of ^+, 
since the non-Gaussian calibr ation to the cosmic var iance 
was derived for this quantity ( Semboloni et al. 20071 ). We 
divide the Poisson term of C++ by two, which compen- 
sates for the additional information of ^_ . The other terms 
contributing to the total covariance (mixed, Gaussian and 
non-Gaussian cosmic variance) do not depend on the num- 
ber of galaxy pairs per bin. Therefore, they are not affected 
by not taking into account the Poisson-noise contribution 
from ^_ and thus they are unchanged. 



6.2. Parameter estimation 

The theoretical model that we fit to the data is a flat 
ACDM cosmology with scale-free, adiabatic and Gaussian 
primordial pert urbations. The t ransfer function is the 
'shape fit' from lEisenstein fc Hu which takes into 

account baryonic suppression; we use a fixed ilb = 0.044. 
The non-linear evolution of the power spectrum is approx- 
imated with the fitting formula of SOS. 

The assumption of scale-invariance is not crucial for 
our results. Indeed, marginalization over the primordial 
spectral index stretches the confidence regions mainly 
along the flm-o's degeneracy direction. The obtained nor- 
malisation for a given Q,^ or the rJm-f 8 relation remains 
unchanged. 

We calculate the log-likelihood on a grid of 6- 
dimensional parameter space: three cosmological parame- 
ters {ilj^,as,h) and three parameters of the redshift dis- 
tribution (a, b, c). 

The Gaussian lensing log-likelihood is 



id, 



(15) 



where an element di of the data vector is either one of the 
measured ^E{0^), {h\^)E{0^) or (M2p)(6'0, and C is the 
covariance of the corresponding estimator. The model nii 
is the theoretical prediction of the shear statistic for the 
same angular separation 0i, and is a function of cosmolog- 
ical and redshift parameters. 

The grid intervals are [0.1; 1] for ilm, [0.4; 1.4] for as 
and [0.6; 0.8] for the Hubble parameter h. The redshift 
parameters values are taken inside of their 2a range: 
[0.5S; 0.69] for a, [6.90; 10.2] for b and [0.49; 0.77] for c. 
Translated into extreme (z) values, this corresponds to 
an exploration range of [0.71; 1.02]. Since the three red- 
shift parameters are correlated, the grid includes models 
that should be rejected by the redshift likelihood alone. 
For this reason we multiply the likelihood, Eq. (fTSj) . by a 
prior given by the likelihood of the redshift distribution 
estimation, 

2 ^ 1 ("■» - n{z,)f 
2 ^ a? 



(16) 



Here rii is the (normalised) number of galaxies in the i- 
th redshift bin of Fig. [5] and n(zi) the fitting function 
Eq. ([M]) . evaluated at the redshift bin centre. The error 
on rii is ai the error bar of the histogram, we neglect the 
cross-correlation between different bins. 

6.3. Constraints 

The left panel of Fig. [TO] shows the marginalised 2D- 
likelihood contours for fl^ and erg using the n{z) of 
Table [TJ A fit to the degeneracy direction yields 

cr8(f7„i/0.25)°'*'5 0.784 ±0.049 for ^e; 

erg (f^m/0. 25)°-^3 = 0.795 ± 0.042 for (|7|^)e; 

(T8(17,n/0.25)"-^'* = 0.785 ± 0.04S for (Mfp). 
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Fig. 10. Left panel: Likelihood contours (1, 2(t) for ilm and as, from the shear correlation function between 1 and 230 
arc minutes (red solid lines), shear top-hat variance between 1 and 230 arc minutes (blue dashed), and aperture-mass 
dispersion between 2 and 230 arc minutes (green dotted-dashed) . A flat, scale-free ACDM model with fib = 0.044 is 
assumed. We marginalise over h and the redshift parameters. Right panel: \a likelihood contours for and cts, from 
the aperture-mass variance between 2 and 35 arc minutes (red solid lines), for scales larger than 35 arc minutes (blue 
dashed) and for scales larger than 85 arc minutes (green dotted-dashed). 



The results for all three statistics are in excellent agree- 



ment . Because of the E-/B-mode mixing ([Kilbinger et al 



200(Th 



we do not use (Mfp) on scales smaller than 
2 arc minutes, therefore we omit the first four data 
poi nts for this statist i c (se e Tabic IB.mB.2l) . If we use 
the Peacock fc Doddi ( 19961 hereafter PD96) non- linear 
power spectrum, the resulting erg is about 2% larger than 
for S03 for a fixed 0.25. 



6.4. Separating small and large scales 

Because of the large connected area of the CFHTLS Wide, 
we are able to obtain interesting cosmological constraints 
by using large scales only. Although the error bars increase 
when small scales are not taken into account, the sensitiv- 
ity to several systematic effects is strongly reduced. The 
deviation from the linear prediction of the shear top-hat 
dispersion is 20% at a scale of 35', for the redshift range 
probed by the Wide survey. The non-linear to linear ratio 
of (M^p) is 3 at 35' and 1.5 at 85', respectively. Our signal 
on large scales is therefore in the linear regime and the 
resulting constraints do not depend on the details of the 
non-linear modeling. In particular, we are not sensitive 
to the difference between PD96 and 803 as can be seen in 
Fig- [HI Other systematics which migh t bias the results on 
small scales are baryonic effects (e.g. Zhan fc Kn ox 2004), 
intrinsic alignment and, maybe most important, shear- 
shape correlations. All these effects are not yet well un- 
derstood as they depend on structure formation on small 
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Fig. 11. The best-fit erg as function of the minimum angu- 
lar scale 0,„in that is used for the x^-analysis. Results are 
shown for S03 (open symbols and Icr error bars) and PD96 
(filled symbols) . The difference between the two non-linear 
models decreases at large scales. 

scales and the relationship between galaxies and dark mat- 
ter. In particular, the shear-shape correlation leads to an 
unde restimation of as ( Hirata fc Seliak|[2004 : Hirata et al. 
20071 ). On scales larger than about 10' the shear field is 



Gaussian. The non-Gaussian calibration of the covariance 
matrix is not needed and also the Gaussian assumption of 
the likelihood is justified. These two factors will yield more 
accurate error estimates on the cosmological parameters. 
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Fig. 12. Left panel: Comparison (l,2,3cr) between our results (bold lines) and the 100 square degree survey (B07, 
filled contours), using in both cases. The redshift distribution is fitted i n the range of [0.2; 1 .5] to be consistent with 
B07. Right panel: Comparison (l,2cr) between WMAP3 (green contours. ISpergel et alJ 2007h and our (Mfp)-results 
between 2 and 230 arc minutes (purple). The combined contours of WMAP3 and CFHTLS Wide arc shown in orange. 



In the right panel of Fig. [10] the results for small and 
large scales are shown. By using only small scales we ob- 
tain tighter constraints than by using only large scales, as 
the signal-to-noise ratio is higher. Using the aperture-mass 
dispersion, the constraints derived from the three angular 
ranges are in very good agreement, with all mean values 
within la: 

as{n,j0.25)°-^^ = 0.780 ± 0.044 for 2' < < 35'; 
c78(f^m/0.25)°-5'^ = 0.780 ± 0.060 for 35' < 61 < 230'; 
cr8(fJm/0.25)"-53 = 0.837 ± 0.084 for 85' < < 230'. 

These results are stable to changes in the smallest an- 
gular scale used. For example, as changes by half a percent 
when only scales larger than 4 arc minutes are used. 

We checked that these constraints are not sensitive to 
possible systematics on angular scales between 50 and 130 
arc minutes, where the B-mode shows a significant bump. 
We fit cosmological parameters using scales with 2' < < 
50' plus 130' < 6 < 230', and found the same results for 
J7m and as- On the other hand, fitting only the affected 
scales, 50' < 6* < 130', we get as = 0.840 ± 0.063 for 
J7rn = 0.25, which is consistent with the results from other 
scales. 




Fig. 13. Comparison (l,2 g) between WMAP3 (green 
contours. ISpergel et al.ll200"7h and our (Mfp)-results in lin- 
ear scale only (85'-230', purple). The combined contours 
of WMAP3 and CFHTLS Wide are shown in orange. 



6.5. Comparison with other data sets 

Our results on cosmological parameters are in very good 
agreement with the most recent cosmic shear analysis 
which combined the first CFHTLS Wide data release, 



the RCS, the VIRMOS-Descart and the GaBoDS sur- 
veys (the '100 square degree survey', B07). In order 
to compare the two results we construct a new Wide 
n(z) histogram that has a consistent redshift distribution. 
Following B07, we only use CFHTLS Deep galaxies with 
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Two-point function 


Angu 


lar scales 






(ItHe 


(l'< 
(2' < 
(2' < 


e < 230') 
e < 230') 
e < 230') 


0.243 ± 0.020 
0.249 ± 0.019 
0.248 ±0.019 


0.771 ±0.030 
0.776 ± 0.029 
0.771 ± 0.029 




(85' < 


e < 230') 


0.255 ±0.027 


0.782 ± 0.038 



a photometric redshift maximum peak probability in the 
ra nge [0.2; 1.5]. We fit an expo nenti al function proposed 



bv iBaugh fc Efstathioul (|l993f ) and IVan Waerbeke et al.. 
(l2002h in the same z-range. The mean redshift (z) — 0.792 



matches the one in B07. The left panel of Fig. [T^ shows an 
excellent agreement between the two results. The compar- 
atively smaller sky coverage of our survey is compensated 
by its larger range of angular scales. It is also interesting 
to notice that our results are in excellent a greement with 
the CTIO survey (jjarvis et al.ll2003l . |2006[) . 

Next, we compare our results for fl^ and erg with 
the Wilkinson Mic rowave Anisotropy Probe 3-year con- 
straints (WMAP3. ISpergel et al.ll2007t ). We combine our 
likelihood with a CMB one computed for a flat ACDM 
cosmology using WMAP3 data only including tempera- 
ture (TT), temperature-polarisation (TE) and polarisa- 
tion (EE) modes. The combination of the two data sets 
leads to remarkably smaller confidence levels as compared 
to individual ones. In particular, as shown in the right 
panel of Fig. [121 the combination of CFHTLS using the 
aperture-mass variance and WMAP3 breaks the severe 
rJni-crg degeneracy. This translates into a reduction of the 
region allowed with 95% confidence level by a factor of 
3.15 as compared to WMAP3 only. The marginalised con- 
straints for each parameter are shown in Table [2l This cor- 
responds to a relative accuracy of 8% in flm and 4% in cts 



improving the WMAP3 constraints of lSpergel et al.l ( 20071 . 
Table 5) by a factor of 1.82 and 1.77 respectively. The 
combinations of CFHTLS and WMAP3 using the shear 
correlation function and top-hat shear variance show con- 
sistent results for and erg as listed in Table [2l 

In view of the weak lensing signal we found on large 
scales, we combine the WMAP3 data with the CFHTLS 
beyond one degree only, and examine the cosmological 
constraints derived from the linear regime. Wc look at 
the constraints on and erg by separating the large an- 
gular scales (85'-230') from the whole sample, which is 
listed in Table O They are shown in Fig. [T21 One can 
see that the large angular scales alone have a significant 
contribution to the total constraint, although the survey 
only covers 57 deg^. It is then realistic to predict from 
this figure that weak lensing surveys may soon be able 
to explore cosmological models using linear theory only, 
similar to CMB physics of primary anisotropics. This is 
very promising for future surveys with sky coverage much 
larger than CFHTLS Wide at the same depth. Equivalent 
constraints from the linear structures, similar to the ones 



shown in Fig. [13] will then be narrower by a factor of at 
least 10. 

Our joint analysis with WMAP3 data is in full agree 
mcnt with similar studies presented in ISpe rgel et al 



(2007,) , using several other data sets. Our estimate for the 
matter dens i ty als o coincides with the result derived by 
lAstier et al. ( 20061 ) based on their SNIa light curves only, 
for a flat ACDM Universe. The comparison with clusters 
of galaxies is, in contrast, less conclusive. Cluster observa- 
tions estimate a broad range of erg values, wi th some be- 
ing fully consistent with our r e sults Gladders et al. (20071). 
(see also iHetterscheidt et al. (2007) for a compilation of 
results), while a recent analysis of simul ations argue for 
higher values ijEvrard et al]l2007t lYepes et al. 2007,) . The 
trends for a high value of erg are a lso derived from a naly- 
ses of the Lyman-alpha forest (see Slosar et al. 20071 . and 
reference therein). 

7. Contamination by shear-shape correlation 

The gravitational lensing signal may be contaminated by 
the intrinsic alignment and by the gravitational shear and 
intrinsic ellipticity (or shear-shape) correlations. We do 
not consider the first term since it would be negligible 
due to a bro ad redshift distribut i on of o ur sa mple. On the 
other hand, Mandelbaum et al. ( 20061 ) and Hirata et al 



(|2007l ) pointed out that the shear-shape anti-correlation 
may bias the estimate of erg by 1 to 20% for a (z) = 1 sur- 
vey on angular scales that we have explored in this work. 
It is therefore important to estimate its amplitude and 
to which extent it may spoil our cosmological constraints. 
We 

We attempt a rather simple analysis of the shear-shape 
correlation (GI) contribution to the shear signal. We use 
the following simple model for the GI correlation func- 
tion £gii which is mo tivated by numerical simulations 
(|Hevinans et al.ll2006bl ) 

(17) 

The lensing efficiency £ is weighted by the source redshift 
distribution 

Xlim Xlim 

5 = y dxin(xi) J dXs'^(Xs) J^^{^) ' 

XI 

For our fiducial flat model with il^ = 0.25 and the redshift 
distribution of Table [H we obtain £ = 95.54 Mpc/h. We 
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fix the scale 6q to 1 arc minute, and further set ^gi = 
on scales larger than 1 degree. 

We perform a combined likelihood analysis using the 
six cosmological parameters as described in Sect. 16.21 and 
the GI amplitude A. The sum ^e+Cgi is fitted to the data. 
Since the 7D-likelihood analysis is very time-consuming, 
we use the marginalised 2a likelihood-region from the pure 
lensing analysis fSect. as a flat prior and do not con- 
sider models outside this region. The marginalised result 
on A is consistent with zero. We find for the amplitude A 
in units of [10^^/i/Mpc arcmin]. 



A = 2.2tli 



for 1' < 61 < 230', 



phological analysis of iMandelbaum et al. I (|2006l) . We do 



p-i 



where the error indicates the 68% confidence region. 
Figure [M] shows there is no significant signal detected at 
any scales. The positive (negative) limit from all scales 
imply a -1-32% (-13%) contamination of the total signal 
by GI at one arc minute. 

Although the confidence region for the constrained GI 
amplitude is large it favours positive correlations, whereas 
from theory we would expect the GI signal to be negative 
( Hirata fc Seliak 2004 ). As a consistency check we used a 
cosmology prior given by the marginalised la likelihood 
region from a pure lensing analysis of the large scale re- 
sults with 6* > 60 arc minutes. The model + Cgi is then 
fitted on scales with 9 < 60 arc minutes. The resulting 
marginalised likelihood for A favours negative GI models 
but is still consistent with zero. This ansatz gives a high 
weight to the large-scale cosmic shear signal, and any sys- 
tematics still present will influence the result. The large 
scale increase in the measured star-galaxy cross correla- 
tion shown in Fig. [7| highlights this concern. As we can- 
not currently distinguish between GI and other possible 
systematic effects we can only conclude from our simple 
analysis that we find no evidence for a non-zero GI signal. 

If our galaxy sample is strongly dominated by high- 
redshift spiral galaxies, then the GI signal may be con- 
siderably weakened, a s one can anticipate from the mor- 



not have enough colour data to explore in detail the spec- 
tral /mor phological types of t he galaxies used in this work. 
However, Zucca et al. ( 20061 ) pointed out that about 80% 
of the VVDS spectroscopic galaxy sample up to i'j^g = 24 
is composed of spiral-like galaxies. It is then possible that 
the fraction of spirals is much higher than elliptical galax- 
ies in the population we are sampling with cosmic shear. If 
so, it would reduc e the contamination to a very small effec- 
tive contribution (jHevmans et al ] l2006bl ). A more detailed 
investigation of the shear-shape analysis using photomet- 
ric redshifts and spectrophotometric information of galax- 
ies is therefore needed and will be discussed in a forthcom- 
ing paper. 



8. Summary and conclusions 

We have presented the weak lensing analysis of the 
CFHTLS T0003 Wide data. The survey covers 57 deg^, 




10 

9 [arcmin] 



100 



Fig. 14. The measured and i^b (open symbols and 
error bars) with the lensing-only best-fit curve (solid blue 
line) and the allowed fractional zblcr-contribution of ^gi 
to the total signal (shaded cyan region). 



about two tirn e s the size of the previous analysis by 
Hoekstra et al. ( 20061 ). and includes a new independent 



field W2. 

The galaxy shape measurements of a homogeneous 
sample of two million galaxies down to = 24.5 have 
been validated using the STEPl and STEP2 simulations 
(jHevmans et all l2006at iMassev et al] l2007bl ). The top- 
hat shear variance, aperture-mass dispersion and the two- 
point shear correlation functions show a significant signal, 
with no galaxy-star correlations, from 1 arc minute up to 
4 degrees. The B-mode is consistent with zero on most 
of these angular scales. It shows, however, a statistically 
significant feature in the range 50-130 arc minutes, of un- 
known origin. We have verified that this feature does not 
influence the cosmological results. 

The two-point statistics show all expected properties of 
a cosmic shear signal up to angular scales 10 times larger 
than the largest non-linear scales of the survey. Hence, 
for the first time the cosmic shear signal can be explored 
with enough confidence to physical scales of about 85 Mpc 
assuming lenses at z = 0.5, for a flat Universe with h = 
0.72 and = 0.27. This is by far the widest scale ever 
probed by weak lensing at that depth. 

The weak le nsing Wide d ata and the photometric red- 
shifts sample of lllbert et al . (2006) are both part of the 
CFHTLS T0003 release and cover common fields. The 
redshift distribution of the Wide data can therefore be 
calibrated using these photometric redshifts, assuming 
with a high confidence level that the two galaxy popu- 
lations are similar. Taking into account the selection cri- 
teria of the weak lensing sample, we find a mean redshift of 
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(z) = 0.949 and a description of the rcdshift distribution 
in excellent agreement with B07. 

Using this redshift distribution, an exploration of con- 
straints on r^m-cTg has been carried out within the angular 
range 1' < < 230'. The marginalised result on ilm-cs 
derived from the aperture-mass variance 

CT8(an/0.25)"-^'* = 0.785 ± 0.043 , 

is in excellent agreement with those obtained by the two 
other statistics (see Fig. [TU| . These constraints perfectly 
match those from B07 as shown in the left panel of Fig. [T^ 
This is interesting because the two samples result from 
complementary approaches: the B07 sample explores the 
consistency of weak lensing results obtained from a het- 
erogeneous sample consisting of four surveys. In contrast, 
our work analyses a very homogeneous data set consist- 
ing of one single survey, and using photometric redshifts 
derived from the same CFHTLS release within the same 
fields as the CFHTLS Wide. 

There is a clear t r end to wards a lower as as conipare d 



to ISemboloni et all (|2006l ) and iHoekstra et"al] (|2006l ) 



This is a result of the less accurate redshift distributions 
used in these analyses which were estimated from the 
Hubble Deep Field photometric redshift sample. This is 
well confirmed when w e combine our pred ictions on 
and as with WMAP3 of lSpergel et al.l (|2007t ). shown in the 
right panel of Fig. [T^l There is a striking diffe rence with 
respec t to the early comparison done by Spergel et al. 
( 2007[) . using the CFHTLS TOOOl results. The 1.5a ten- 
sion is no longer visible; in contrast, there is a large overlap 
between the two data sets. The joint CFHTLS- WMAP3 
likelihood analysis then leads to tight marginalised con- 
straints on and erg, 



0.248 ± 0.019 and erg = 0.771 ± 0.029 



corresponding to an accuracy of 8% and 4% on these two 
parameters. Hence, using a much better photometric red- 
shift sample, based on the Deep CFHTLS T0003 data sets 
that directly calibrate the genuine CFHTLS galaxy popu- 
lation, removes one of the primary uncertainties of earlier 
CFHTLS weak lensing analysis. 

Considering the potential nuisances of systematic ef- 
fects related to non-linear scales, we split the sample 
into three ranges of angular scales: the 'highly non-linear' 
(2' < < 35'), the 'intermediate' (35' < 9 < 230') and 
the 'linear' (85' < < 230') scales. The analysis of the 
three sub-samples do not reveal significant differences be- 
tween each regime (see Fig. [101 right panel). The results 
are also stable to changes in the lower angular scales in- 
creasing from 2' to 4'. This shows that the CFHTLS Wide 
cosmic shear survey is not yet dominated by uncertainties 
related to our poor knowledge of astrophysical systematics 
at small scales. Finally, we find that excluding scales with 
a significant B-mode (50' < < 130') from the analysis 
does not change our results. The constraints on flm — cr^ 
are therefore insensitive to the level of residual systemat- 
ics in our data. All these tests strengthen the confidence 
and reliability of our results. 



The very large range of angular scales explored by the 
CFHTLS Wide opens a new window to cosmic shear sur- 
veys. It enables for the first time a comparison of cosmic 
shear and WMAP3 signals using only linear scales. The 
constraints shown in Fig. 1 131 demonstrate that there is still 
great predictive power from the linear regime only. Future 
weak lensing surveys which cover areas significantly larger 
than the CFHTLS will be able to pin down a much nar- 
rower region in parameter space. Thus, it will be possible 
to obtain cosmological parameters to percent-level accu- 
racy and below from combining CMB and weak lensing 
using linear theory. 

Finally, the impact of the contamination by the shear- 
shape correlation on cosmic shear surveys like CFHTLS 
is still unclear. We find its amplitude to be very low and 
compatible with zero a t all scales w e explo red. The low 
amplitude derived from iHirata et al.l (2007), using a sur- 
vey shallower than the CFHTLS Wide, had already sug- 
gested that it should be a small effect and a difhcult-to- 
detect signal in the CFHTLS- Wide, in particular if our 
gal axy sample is dom inated by high-redshift spiral galax- 
ies ( Zucca et al. 20061 ). At present, we can measure erg to 
a precision of about 5% and so this bias is still reasonably 
low. With future work, however, this bias may become the 
main source of error. 

The CFHTLS is still in progress and the next re- 
lease will include more sky coverage and also a new field, 
W4. In this work, we only use the wide z'-band data to- 
gether with the photometric redshift from the Deep T0003 
u* , g' ,r' ,i' , z' . The next step is therefore a more detailed 
analysis of multi-colour data sets. A better check of sys- 
tematics will be possible by cross-correlating the lensing 
signal obtained independently in indifferent filters. The 
larger CFHTLS Wide sample with u*,g', r',i', z' wiU also 
improve tomographic studies and will provide a photo- 
metric redshift to each individual galaxy. We will then 
be in the position to better control contaminations by in- 
trinsic alignment an d the shear-shape (GI) correlations 
( Bridle fc Kind 2007) and to move towards a full tomo- 
graphic exploration of the CFHTLS Deep and Wide sur- 
veys together. 
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Appendix A: STEP simulation calibration 

One of the crucial issues for weak lensing studies is the re- 
liability of galaxy shape measurement and the control of 
systematics. The detection and measurement of weak lens- 
ing is a technical challenge. Weak distortion induced by 
gravitational lensing in the observed shapes of galaxy im- 
ages is only ~ 1%, much smaller than the typical intrinsic 
ellipticity dispersion ~ 30%. To further comphcate the sit- 
uation the observed shape of the galaxies i s affected by the 
PSF. The Shear TEsting P rogramme^ (Heyma ns et al 



2006al: iMassev et all l2007bl) . hereafter STEP, is a col- 



laborative project aiming to calibrate and improve weak 
lensing methods using realistic Wide field simulated im- 
ages. The first and second generation of STEP simula- 
tions (hereafter STEPl and STEP2) are designed for a 
ground-based survey. In order to check the reliability of 
the shear measurement used in this analysis, we calibrated 
the pipehne using all data sets from STEPl and STEP2. 

STEPl simulations contain relatively simple galaxy 
morphologies generated using the SkyMaker software 
Five constant shears, 7^"^"* =[0.0, 0.005, 0.01, 0.05, 0.1], 
are applied to the galaxies, while the second component 
^mput always set to zero. Finally, galaxy and stellar 
point sources are convolved with six different constant 
PSFs which attempt to reproduce PSF shapes, that are 
typical of ground-based observations. In this way 30 sets 
of images, differing in PSF type and/or shear strength are 
produced. Each set is composed of 64 images. The sky 
noise is spatially uncorrelated. 

STEP2 simulations contain complex galaxy mor- 
phologies produced using a shapelet simulation package 



(jMassev et al.l 120041 ). Six sets of 64 images with random 



constant input shears are convolved each with a different 
optical PSF. The six PSFs are chosen to span a range 
of realistic ground-based observing conditions. For each 
image, a twin image is produced, in which galaxies are ro- 
tated by 90° before applying the same shear and the same 
PSF. Combining the shear analysis on rotated and non- 
rotated images demonstrates the pure measurement bias, 
since the noise due to the scatter in a galaxies' intrinsic 
morphology is removed. The model of the sky noise is also 
more complex than the one adopted to generate STEPl 
simulations. It is in fact a correlated noise which aims to 
reproduce the noise of the drizzling process. 

Our pipeline is an application of the KSB-I- method. 
The observed galaxy shape is modeled as a convolution of 
the sheared galaxy with the PSF, which in turn is modeled 
as a circular profile convolved with a small anisotropy. 
Assuming the mean of the intrinsic ellipticity distribution 
of galaxies to be zero and the PSF anisotropy to be small, 
the first-order of the shear, 7, can be computed from the 
observed ellipticities of galaxies, e"^** as follows: 



(A.l) 



9 
10 



http : //www .physics .ubc . ca/heymans/ step . html | 
http : / / terapix . iap . f r/ soft/ skymaker 
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IS 



where P**™ is the smear polarisability and q — 
the anisotropic component of the PSF. The symbol * in- 
dicates those quantities a. r e mea sured on stars. Pj, de- 



fined in lLuppino fc Kaisen (|l997t ). is the correction to the 



shear polarisabihty which includes circular smearing by 
the PSF. 

We compute the stellar quantities, p^^* ^ psm* g^j^^ 
with the same filter function W {9, a) in o rder to keep the 



calibr ation free of extra bias. F ollowing iHoek stra et al. 
(Il998h and the ST EP results (|Hevmans et ah ,2006a; 
Massev et al. 2007bl) . for each galaxy we compute all 



quantities, including those estimated from the stars, using 



a filter scale a 



as given by IMC AT. 



We did not apply the same PSF anisotropy correction 
to small and large objects. Using the STEPl and STEP2 
simulated catalogues we found that the measurement of 
moments from small objects can be significantly improved 
and are more robust by first resampling the intensity of 
light in each pixel. Each image is oversampled by a factor 
of two and interpolated using a nearest neighbour inter- 
polation kernel, prior to measure shapes of objects. The 
interpolation works very well for objects with a size close 
to the star size and does not produce any detectable extra 
bias. However, it fails and may even degrade the signal 
as the object size increases. STEP simulations show the 
transition arises when object size exceeds 1.2 x seeing. 

We approximate by half of its trace, TrP^/2. Since 

individual TrPy are noisy, we derive their values from a 

fit as function of some 
r 



in Hev mans et al.l ( 2006. 



galax y properties. As desc ribed 
^ and lMassev et all (|2007bh . the 



shear bias parameter m often depends on object sizes, r^, 
and magnitudes mag. We therefore fit TrP^ in the rg- 
mag plane using a polynomial that only depends on these 
two parameters. The Tr P-y-dependence on mag is more 
scattered than that on rg, so we choose a function that 
gives more weight on r^: 



multiplicative calibration bias mi for each of the STEPl 
PSF models of the simulation. Averaging over the six 
STEPl PSF models, our shear measurement bias is less 
than 1% as can be seen in Fig. IA.2I (left panel). The two 
right panels of Fig. I A. 2 1 show the bias found in the STEP2 
simulations, once the rotated and unrotated image s have 
been merged as described in iMassev et all (|2007b[ ). The 
two components of the shear are underestimated by about 
3% on average. Model C is the most similar to the seeing 
found in our CFHTLS images. It is worth noting that our 
poorest results come from PSF 2 of STEPl and PSFs D 
and E of STEP2 which have the strongest anisotropy of 
all the simulations. Every shear method tested on these 
particular simulations had difficulty recovering the correct 
shear. 

These results show the residual bias is well constrained 
and reasonably low for our purposes. In fact, the shape 
measurement bias is much lower than the total error af- 
fecting the cosmological parameter estimation. It should 
also be noticed that the STEPl simulations have charac- 
teristics much more similar to the Wide data than STEP2 
simulations. In particular, the structure of the noise and 
the PSF types adopted to generate the STEPl simulations 
are very close to the ones of our data. For this reason the 
effective bias in the CFHTLS Wide shear catalogue is ex- 
pected to be about 1%. 

Appendix B: Shear two-point correlation data 

The data vectors and Icr error bars plotted in Fig. |4l 
for the various shear two-point functions are listed in 
Table IrT1IR2] 



TrP 



= ai + a2 Tg + as Tg + a^ mag 



(A.2) 



Table ETT] summarises the key parameters of our shear 
measurement pipeline. As an illustration of its application, 
the left panel of Fig. lA.ll shows a compilation of all PSF 
measurements for all stars of the 57 pointings used in this 
work. The distribution of corrected stars ellipticities in the 
right panel, shows a reduction by a factor of 10 in both 
the average ellipticity and dispersion, without showing a 
preferential direction. 

Applying the pipeline to the STEP simulations, we 
quantify the STEP result s using th e fit defined in 



Hevmans et al.l (|2006al ) and lMassev et al.l (j2007bf ). which 
expresses the difference between measured and input shear 
through a linear relation: 



(7.: 



input input , 



(A.3) 



where « = 1, 2 are the two shear components. For a perfect 
shear measurement, and Ci would be zero. Figure [A. 2! 
shows the values of the residual shear oflFset ci and of the 
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Source Detection 


hfindpeaks 




PSF: 2D polynomial model 


2nd order fit of e*{rg), P"'"*(rg) and P"''*(rg) 




Galaxy radius size 


Vg from hfindpeaks 




Quadrupole, P"" and P"^ estimate 


Interpolation (seeing < Vg < 1.2 x seeing ), 6'max = ^rg 


and AO = 0.5 pixel 




Approx {rg> 1.2 x seeing ), Omax ~ Int[4r3] and AO 


— 1 pixel 


correction 


Fit in {rg,mag) to Tr(P^)/2 




Weight 


Hoekstra et al. f2002bl 




EUipticity cut 
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Size cut 
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Significance cut 
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Fig. A.l. Left panel: The observed ellipticities of all stars in the 57 pointings. The red cross marks the (0, 0) position. 
Right panel: The residual star ellipticity after PSF anisotropy correction. In both plots, the mean values of the two 
ellipticity components are given. 
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Fig. A. 2. The calibration bias m and the residual offset c of our pipeline estimated using STEP simulations. Left 
panel: The results of STEPl for the first component of shear. PSF models are labeled from to 5. Middle and right 
panels: The results of STEP2 for the two shear components. PSF models are labeled from A to F. 
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Table B.l. Values of the shear correlation function and the shear top- hat variance, as function of scale 9 in arcmin. 
The errors include statistical errors and non-Gaussian-calibrated cosmic variance for the E-mode, while only statistical 
uncertainty contributes to the error of the B-mode. 



e 




& 








<5Ce 






5Cb 






(|7P>B 


5{b?U 




1.00 


9 


704e- 


05 


2.494e-05 


3 


416c- 


05 


3 


200C-05 


1 


158C-04 


-1.809C-06 


1 


749c- 


05 


1 


036C-05 


1.18 


8 


548e- 


05 


2.546e-05 


2 


947c- 


05 


2 


739C-05 


1 


054e-04 


-2.083C-06 


1 


566c- 


05 


8 


874C-06 


1.39 


6 


511e- 


05 


7.956e-06 


2 


555e- 


05 


2 


352e-05 


9 


631e-05 


-1.422e-06 


1 


405e- 


05 


7 


606e-06 


1.64 


6 


160e- 


05 


-1.138e-05 


2 


209e- 


05 


2 


012e-05 


8 


777e-05 


-1.570e-07 


1 


261e- 


05 


6 


521e-06 


1.93 


7 


294c- 


05 


-1.479C-05 


1 


925c- 


05 


1 


730C-05 


7 


933C-05 


9.634C-07 


1 


133c- 


05 


5 


593C-06 


2.28 


6 


865e- 


05 


-2.900e-06 


1 


672e- 


05 


1 


486e-05 


7 


llle-05 


1.695e-06 


1 


020e- 


05 


4 


800e-06 


2.69 


4 


578e- 


05 


1.275e-05 


1 


456e- 


05 


1 


277e-05 


6 


373e-05 


2.020e-06 


9 


174e- 


06 


4 


122e-06 


3.17 


4 


167e- 


05 


6.394e-06 


1 


279e- 


05 


1 


lOle-05 


5 


689e-05 


2.544e-06 


8 


276e- 


06 


3 


542e-06 


3.74 


2 


826e- 


05 


2.609e-06 


1 


117e- 


05 


9 


473e-06 


5 


lOOe-05 


2.656e-06 


7 


459e- 


06 


3 


045e-06 


4.41 


4 


490c- 


05 


-1.043C-06 


9 


840c- 


06 


8 


201C-06 


4 


615C-05 


2.538e-06 


6 


735c- 


06 


2 


617e-06 


5.20 


2 


325e- 


05 


8.105e-06 


8 


646e- 


06 


7 


050e-06 


4 


190e-05 


2.210e-06 


6 


081e- 


06 


2 


247e-06 


6.13 


2 


908e- 


05 


-9.270e-07 


7 


556e- 


06 


5 


996e-06 


3 


756e-05 


1.671e-06 


5 


487e- 


06 


1 


929e-06 


7.22 


2 


979e- 


05 


1.723e-06 


6 


634e- 


06 


5 


099e-06 


3 


346e-05 


1.241e-06 


4 


962e- 


06 


1 


656e-06 


8.52 


3 


179e- 


05 


-5.419e-06 


5 


846e- 


06 


4 


398e-06 


2 


986e-05 


1.020e-06 


4 


483e- 


06 


1 


421e-06 


10.04 


1 


668c- 


05 


-4.458C-06 


5 


195c- 


06 


3 


793C-06 


2 


666e-05 


7.947C-07 


4 


070c- 


06 


1 


219C-06 


11.84 


1 


687e- 


05 


2.365e-06 


4 


647e- 


06 


3 


257e-06 


2 


393e-05 


5.714e-07 


3 


711e- 


06 


1 


044e-06 


13.97 


1 


688e- 


05 


1.362e-06 


4 


189e- 


06 


2 


748e-06 


2 


154e-05 


4.431e-07 


3 


408e- 


06 


8 


937e-07 


16.47 


1 


530e- 


05 


-1.318e-06 


3 


830c- 


06 


2 


338e-06 


1 


944e-05 


2.889C-07 


3 


159c- 


06 


7 


649e-07 


19.42 


1 


579e- 


05 


-3.319C-06 


3 


530c- 


06 


2 


004e-06 


1 


751e-05 


1.432C-07 


2 


942c- 


06 


6 


548G-07 


22.90 


1 


353c- 


05 


7.628C-07 


3 


296e- 


06 


1 


727e-06 


1 


578e-05 


-8.315e-08 


2 


763c- 


06 


5 


609C-07 


27.00 


1 


207e- 


05 


-5.854e-07 


3 


051e- 


06 


1 


477e-06 


1 


413e-05 


-3.782e-07 


2 


600e- 


06 


4 


812e-07 


31.84 


9 


731e- 


06 


-1.522e-06 


2 


850e- 


06 


1 


274e-06 


1 


266e-05 


-7.709e-07 


2 


453e- 


06 


4 


135e-07 


37.54 


1 


057e- 


05 


-2.188e-06 


2 


696c- 


06 


1 


llOc-06 


1 


138C-05 


-1.175C-06 


2 


316c- 


06 


3 


562C-07 


44.26 


6 


947c- 


06 


-2.544e-06 


2 


513c- 


06 


9 


671C-07 


1 


007C-05 


-1.595C-06 


2 


179c- 


06 


3 


076C-07 


52.19 


9 


153c- 


06 


-5.789C-06 


2 


363c- 


06 


8 


417e-07 


8 


614C-06 


-1.932C-06 


2 


051e- 


06 


2 


664C-07 


61.54 


7 


506e- 


06 


-3.359e-06 


2 


239e- 


06 


7 


335e-07 


7 


216e-06 


-2.172e-06 


1 


922e- 


06 


2 


315e-07 


72.57 


4 


613e- 


06 


-4.799e-06 


2 


081e- 


06 


6 


492e-07 


5 


996e-06 


-2.245e-06 


1 


799e- 


06 


2 


020e-07 


85.57 


1 


110c- 


06 


-2.967C-06 


1 


977c- 


06 


5 


706C-07 


4 


980C-06 


-2.150C-06 


1 


680c- 


06 


1 


774C-07 


100.90 


2 


006e- 


06 


-3.665e-06 


1 


873e- 


06 


5 


117e-07 


4 


159e-06 


-1.927e-06 


1 


567e- 


06 


1 


570e-07 


118.98 


2 


416e- 


06 


-2.401e-06 


1 


797e- 


06 


4 


636e-07 


3 


619e-06 


-1.689e-06 


1 


465e- 


06 


1 


407e-07 


140.29 


1 


982e- 


06 


-1.174e-06 


1 


743e- 


06 


4 


272e-07 


3 


203e-06 


-1.482e-06 


1 


373e- 


06 


1 


284e-07 


165.42 


1 


672e- 


06 


-4.094e-07 


1 


727e- 


06 


4 


059e-07 


2 


745e-06 


-1.171e-06 


1 


301e- 


06 


1 


206e-07 


195.06 


3 


439c- 


06 


-9.141C-07 


1 


743c- 


06 


3 


961C-07 


2 


214C-06 


-7.842C-07 


1 


259c- 


06 


1 


223C-07 


230.00 


1 


780e- 


06 


-1.420e-06 


1 


876e- 


06 


4 


081e-07 


1 


615e-06 


-2.919e-07 


1 


333e- 


06 


1 


686e-07 



L. Fu et al.: Very weak lensing in the CFHTLS Wide 
Table B.2. Values of the the aperture-mass variance, as function of scale 9 in arcmin. 



e 






(Mi) 






S{Ml) 


1.00 


1 


201e-05 


1.143e-06 


5 


944e-06 


5 


695e-06 


1.18 


1 


402e-05 


-3.377e-06 


5 


138e-06 


4 


875e-06 


1.39 


1 


448e-05 


-3.000e-06 


4 


447e-06 


4 


173e-06 


1.64 


1 


306C-05 


-8.902C-07 


3 


856e-06 


3 


571C-06 


1.93 


1 


107e-05 


-2.830e-07 


3 


353e-06 


3 


059e-06 


2.28 


9 


292e-06 


-8.651e-07 


2 


923e-06 


2 


622e-06 


2.69 


8 


821e-06 


-1.441C-06 


2 


552C-06 


2 


248C-06 


3.17 


8 


556C-06 


-1.258C-06 


2 


238C-06 


1 


930C-06 


3.74 


8 


472e-06 


-9.641e-07 


1 


973e-06 


1 


657e-06 


4.41 


8 


221e-06 


-9.198e-07 


1 


743e-06 


1 


423e-06 


5.20 


7 


676e-06 


-7.686e-07 


1 


542e-06 


1 


222e-06 


6.13 


6 


249C-06 


-5.371C-07 


1 


364e-06 


1 


048e-06 


7.22 


5 


030C-06 


-1.800C-07 


1 


209e-06 


8 


971e-07 


8.52 


4 


609e-06 


3.606C-07 


1 


079C-06 


7 


693e-07 


10.04 


4 


508e-06 


5.421e-07 


9 


612e-07 


6 


611e-07 


11.84 


4 


301e-06 


4.289e-07 


8 


544e-07 


5 


684e-07 


13.97 


3 


976e-06 


3.853C-07 


7 


688C-07 


4 


874C-07 


16.47 


3 


526C-06 


3.019C-07 


7 


046C-07 


4 


171C-07 


19.42 


3 


009C-06 


2.118C-07 


6 


501C-07 


3 


572C-07 


22.90 


2 


600e-06 


1.249e-07 


5 


954e-07 


3 


066e-07 


27.00 


2 


239e-06 


4.597e-08 


5 


387e-07 


2 


633e-07 


31.84 


1 


976e-06 


2.932C-08 


4 


840C-07 


2 


261C-07 


37.54 


1 


877e-06 


8.679e-08 


4 


359e-07 


1 


944e-07 


44.26 


1 


675e-06 


1.961e-07 


3 


967e-07 


1 


674e-07 


52.19 


1 


322e-06 


3.382e-07 


3 


668e-07 


1 


443e-07 


61.54 


1 


147e-06 


4.524e-07 


3 


454e-07 


1 


245e-07 


72.57 


1 


303e-06 


5.266C-07 


3 


312e-07 


1 


077e-07 


85.57 


1 


475e-06 


5.585e-07 


3 


223e-07 


9 


356e-08 


100.90 


1 


511e-06 


4.750e-07 


3 


172e-07 


8 


174e-08 


118.98 


1 


476C-06 


2.712C-07 


3 


158C-07 


7 


206e-08 


140.29 


1 


317C-06 


3.462C-09 


3 


174C-07 


6 


437C-08 


165.42 


9 


472C-07 


-9.410e-08 


3 


220e-07 


5 


879e-08 


195.06 


5 


901e-07 


-1.055e-07 


3 


345e-07 


5 


692e-08 


230.00 


3 


641e-07 


-7.764e-08 


3 


966e-07 


7 


353e-08 



